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Abstract 

We consider a {d + 2)-dimensional class of Lorentzian geometries holographically dual to a rela- 
tivistic fluid flow in (d+1) dimensions. The fluid is defined on a (<i+l)-dimensional time-like surface 
which is embedded in the (d + 2)-dimensional bulk space-time and equipped with a flat intrinsic 
metric. We find two types of geometries that are solutions to the vacuum Einstein equations: the 
Rindler metric and the Taub plane symmetric vacuum. These correspond to dual perfect fluids with 
vanishing and negative energy densities respectively. While the Rindler geometry is characterized 
by a causal horizon, the Taub geometry has a timelike naked singularity, indicating pathological 
behavior. We construct the Rindler hydrodynamics up to second order in derivatives of the fluid 
variables and show the positivity of its entropy current divergence. 
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I. INTRODUCTION 



The holographic principle proposes that (d + 2)-dimensional (quantum) gravitational 
theories are equivalent to (d + l)-dimensional field theories living on a boundary of the 
higher dimensional spacetime. The most concrete examples of holography are the so-called 
gauge/gravity correspondences, where quantum gravity with negative cosmo logical constant 
is dual to certain flat spacetime gauge theories. The gauge theory can be thought of as living 
on the timelike boundary at spatial infinity, in which the bulk spacetime is holographically 
encoded. A particularly interesting consequence of this duality is that the hydrodynamics 
of the gauge theory can be effectively described by the long time, long wavelength dynamics 
of a black hole living in the bulk. In this fluid-gravity correspondence, Q, [2] the Navier- 
Stokes equations of the fluid are equivalent to the subset of the General Relativity (GR) field 
equations called the momentum constraints, which constrain data on the timelike boundary 
surface. 

However, by studying the dynamics of a Rindler acce!e r arion horizon in flat spacetnne g 
|5j, two of the authors argued that the relationship between holography and hydrodynamics 
is not limited to theories with negative cosmological constant. In the Rindler wedge of 
flat spacetime, the usual Minkowski vacuum is a thermal state at finite temperature (see 
for example [6|). Just as in the AdS/CFT examples, one can study the long wavelength, 



rich are expected to be dual to 
3, |8] it has been shown that one 



long time perturbations of this background spacetime, w 
the hydrodynamics of the thermal state. Indeed, recently 
can construct explicit bulk solutions to the vacuum Einstein equations dual to a particular 
non-relativistic fluid by perturbing around the Rindler geometry. The holographic fluid 
in this case is defined on an arbitrary timelike surface S c of fixed radial coordinate r = 
r c in the bulk geometry {9]. These are the hyperbolas associated with the worldlines of 
accelerated observers. Working in a non-relativistic hydrodynamic expansion, one can solve 
the field equations subject to the boundary conditions of a fixed flat induced metric on S c 
and a regular event horizon. The momentum constraints on S c again are the non-relativistic 
Navier-Stokes equations. 

This result implies that there is an underlying duality between a field theory on S c and 
the bulk interior Rindler space. The nature of holography in asymptotically flat spacetimes 
has remained mysterious and it expected that the dual field theory will be nonlocal 10 ]. 
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An intriguing aspect of these results is that the dual fluid thermodynamics, constructed 
from the quasi-local Brown- York stress tensor, is characterized by zero equilibrium energy 
density even though there is non-zero temperature. Similarly, one can also show the entropy 
density of the fluid, equivalent to the area entropy density of horizon, is independent of its 
temperature. 

Despite these unusual thermodynamical properties, the hydrodynamics of the Rindler 
fluid appears to be well-defined and contains no obvious pathologies. Studying this fluid 
in more detail may yield additional clues into the nature of the microscopic duality. In Q 
the authors showed that their results can be obtained as the non-relativistic limit of an 
underlying relativistic fluid. In particular, they constructed a general theory for a viscous 
relativistic fluid with zero energy density and by matching to the non-relativistic solution 
were able to determine some of the viscous transport coefficients. While shear viscosity 
to entropy density ratio saturates the Kovtun-Son-Starinets bound of l/4ir llj, the bulk 
viscosity is not an independent transport coefficient even though the fluid is non-conformal. 
In the second order viscous hydrodynamics, there are six independent transport coefficients, 
but by matching to the non-relativistic solution one is only able to determine four of these. 



Later, 



12] studied how higher derivative corrections to the gravitational field equations 



affect the properties of the dual fluid. In the AdS/CFT correspondence, such terms are 
associated with quantum corrections or other deformations, which modify the values of the 
transport coefficients 13j. Interestingly, in this case only the second (and higher) order 
hydrodynamics is affected; the shear viscosity to entropy density ratio and the first order 
Navier-Stokes equations are universal. 

In this paper, our main goal is to expand upon the results of js| by completely determining 
the relativistic fluid dual to the Rindler spacetime. We start by considering a particular class 
of (d + 2)-dimensional Lorentzian geometries. These metrics are stationary and on the slices 
S c the intrinsic metric is flat and the extrinsic curvature is such that the Brown- York quasi- 
local stress tensor has a perfect fluid form. Thus, these metrics can, in principle, act as 
the bulk gravitational dual to a fluid on S c . Solving the vacuum Einstein equations, we 
find there are two possible branches of solutions. One is the Rindler solution described 
above and the other is the known as Taub plane symmetric vacuum and is associated with a 
fluid of negative energy density. The Taub metric has non-trivial curvature and, crucially, a 
naked singularity, which indicates pathological behavior in the dual field theory. In contrast, 
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the Rindler solution is well-behaved, and using the formalism developed in [8j we find the 
bulk solution and fluid stress tensor in a derivative expansion up to second order, fixing the 
remaining two transport coefficients. In the fluid-gravity correspondence, the fluid entropy 
current is mapped into the area current [l^] of the evolving horizon surface. We compute 
this current to second order and find that its divergence is positive definite, consistent with 
Hawking's area theorem. 

The plan of this paper is as follows. In Section II, we review the general construction of 
the solutions developed in j?-9| and describe the Rindler and Taub geometries. In Section 
III, we present the earlier non-relativistic description of the Rindler fluid and use this to 
develop and eventually calculate the fully relativistic metric and its corresponding stress 
tensor. Section IV is devoted to the calculation of the holographic area entropy current and 
its divergence. We then conclude with a brief discussion and description of possible future 
work. 



II. RICCI FLAT GEOMETRIES AND FLUIDS 

In the following we will construct certain (d + 2)-dimensional Lorentzian geometries holo- 
graphically dual to a fluid flow in (d + 1) dimensions. The fluid is defined on a (d + 1)- 
dimensional timelike surface S c embedded in the (d + 2)-dimensional bulk space-time. We 
choose the timelike surface to be defined by fixed bulk radial coordinate, r = r c . Consider 
the following metric ansatz for the bulk geometry [si 

ds 2 = g AB dx A dx B = -h{r)dt 2 + 2dtdr + e^Wdx* , (1) 

where x A = (t,x\r), i = l..d and d>2. On surfaces of r = r c , where r c is a constant, the 
induced metric is 

ds 2 = lllv dx»dx v = -h(r c )dt 2 + e 2r{rc) dxidx* . (2) 

This metric is flat, as can be seen by the coordinate re-scaling i = \Jh{r c ) and x l = e T ^ Tc ^x' l } 
which leads to the standard Minkowskian form 

ds 2 = jp i pdx' 1 dx L ' = —dt 2 + dxidx 1 , (3) 

where x^ = (t, x l ) 
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The Brown- York stress-energy tensor [15( (in units where 167rG = 1) associated with the 
r = r c slice is 

= 2{K lilv - K^) , (4) 

where K^ v = ^CnI^u and Cn is the Lie derivative along the normal to the slice N A . Using 
([T]) we find that 

T E = p = -2dVhr', Tj-j = p = 2Vh ({d - l)r' + % , (5) 

where primes represent derivatives with respect to r and the expressions are evaluated at 
r = r c . The stress-energy tensor has the form of a perfect fluid with energy density p and 
pressure p. 

We wish to determine whether there is a solution to the vacuum Einstein equations 

Rab = 0, (6) 

of this general form. The Hamiltonian constraint, RabN a N b , where N A is the unit spacelike 
normal to the r = r c slices, is 

G AB N A N B = TZ- K^K^ + K 2 , (7) 

where 1Z is the Ricci scalar associated with the induced metric of the slice. Re-expressing 
this equation in terms of the Brown- York stress tensor and using the fact that 1Z — for 
r = r c , we get 

dT^T^ = T 2 . (8) 

Inserting the general form of a perfect fluid stress tensor and (jSJ), one finds that this condition 
is satisfied by two types of equations of state [8j 

—2d 

(z) p = 0, (ii) p =—-—p. (9) 

A. The p = case: Rindler Geometry 

Consider first the equation of state p = 0. Using Eqn. (jSJ), this condition implies that 
t' — and as a result, r = const.. The remaining field equations imply that h(r) = r, and 
we get the metric 

ds 2 = -rdt 2 + 2dtdr + dxidx 1 , (10) 



which describes a region of flat (d + 2)-dimensional Minkowski space-time in "ingoing 
Rindler" coordinates. The null surface r = acts as a causal horizon to accelerated ob- 
servers, whose world-lines correspond to the surfaces of constant r = r c . Although the 
Rindler metric is just a patch of flat space-time, the associated quantum field theory on this 
background has many of the same properties as a black hole solution due to the existence 
of the causal horizon. In particular, surfaces of r = r c have a local Unruh temperature (in 
units where H — c — 1) 

T=—!—. (11) 



c 

Strictly speaking, a Rindler horizon does not have a Bekenstein-Hawking entropy density. 
However, one can assign the Rindler horizon this entropy based on the holographic principle, 
or, more concretely, take the entropy to be the thermal entanglement entropy of the quantum 
fields in Rindler wedge [if]] . This statistical entropy scales like an area, but is a UV divergent 
quantity. If a Planck scale cutoff is chosen appropriately, the entanglement entropy agrees 
with the Bekenstein-Hawking formula, i.e. in units where \6ttG = 1 

s = 4tt . (12) 

Given the existence of an equilibrium Unruh temperature and a Bekenstein-Hawking entropy 
density, the metric ( ITUl) may be considered as providing a dual geometrical description of 
a perfect fluid with zero energy density in one lower space dimension. We will discuss the 
hydrodynamics of this case in detail in section III. 

B. The p < case: Taub Geometry 

Consider next the second equation of state in ([9]). In order to construct the background, 
we plug-in the values of the energy density and pressure in terms of the metric functions 
into the equation of state, which gives 

^ih'h - < 13 > 

Consider the equation R rr = 0. It yields 

r' 2 + t" = , (14) 
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which is solved by , 

r(r) = ln(dr + C 2 ) . (15) 

Inserting this into f lT3|) . we find 

k{r) = Wr + C 2 ) d ^ • (16) 
Therefore the dual gravitational solution is 

d s 2 = -— — ° 3 dt 2 + Idtdr + (Cir + C 2 fdx i dx i . (17) 
[C\r + 62) 

Redefining the radial coordinate f = Cir + C2 and re-scaling the time coordinate, this metric 
takes the form 



A 



ds 2 = -—jzidt 2 + Idtdf + f 2 dxidx l , (18) 
where A is a constant. 

I — I 

In four-dimensions (d = 2) this metric was found by A. H. Taub in 1951 [17j. It can be 
considered the vacuum solution exterior to an infinite plane-symmetric object with uniform 
mass density. The Kretschmann scalar for this solution reads 

RabcdR ABCD ~ ~^2^+Tj ' ( 19 ) 

which implies that there is a curvature singularity at f = and the solution is asymptotically 
flat at infinity f = 00. The curvature singularity at f = is timelike and naked, consistent 
with the fact that the energy density computed from the Brown- York stress tensor 



-2dVA 

p = (20) 



is always negative. The global structure of this metric was analyzed in [18| . At infinity there 
are two flat null surfaces, while the timelike naked singularity is located in the interior. 

While it seems clear that this branch is problematic, let us nevertheless make a few 
remarks. 

(i) One can make a spatial boost and re-write the metric in terms of the energy density p. 
This yields 

ds 2 = - P ^--u il u v dx ii dx v + 2uadx^df + f 2 P„ v dx^dx v . (21) 
4cr 
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One could then formally allow the variables p(x^) and u^(x^) and solve the field equations 
order by order in a derivative expansion in d^p and d^Uy as is done in the fluid-gravity 
correspondence jlj]. However, unlike the Rindler solution, in this case there is no casual 
horizon in the background zeroth order solution. Therefore, interpreting this geometry 
as being dual to a finite temperature state of a field theory (and then perturbations of this 
state as hydrodynamics) is problematic. For instance, imposing the thermodynamic identity 
p + P = sT yields 

sT=(t±±)p. (22) 



l-d / 

Since p > 0, this equation implies formally a state of negative temperature (or negative 
entropy). A related fact is that the squared speed of sound v 2 s = ^ is negative, which 
indicates the dual field theory is characterized by an instability. It would be interesting to 
understand the role of this type of exotic solution in asymptotically flat holography, (ii) 
Following 8[ we can create a scalar field Lagrangian that mimics the equation of state for 
the Taub metric. We consider 

1 = J d d+l ^-gF{X,<j>) , (23) 

where X = —(l/2)g fJiU d fJ- <j)d u <j). The stress-tensor is given by 

dF 

T, v = -2— + g, u F . (24) 



This gives 



dF 

T» u = -^d,<pd u <P + g, u F . (25) 



If we identify the four-velocity (of a potential flow) 

a 



•*-m> (26) 

the stress-tensor has the form of a perfect fluid with pressure F and 

P = 2A-§-F. (27) 

Imposing the equation of state, we find the condition on F(X), which leads to an action of 
the form 

I = j d d+lx V^X-^ . (28) 
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III. THE RINDLER/FLUID CORRESPONDENCE 



A. General Setup 

In order to study the hydrodynamics of the fluid living on the r = r c slices in Rindler 
geometry, we have to perturb this background. The metric (TT0|) corresponds to a fluid in 
thermodynamical equilibrium in the rest frame of some observer. The first step is to make a 
set of coordinate transformations to obtain a new metric corresponding to the fluid flow in 
an arbitrary frame. More geometrically, these transformations keep the induced metric at 
r c flat, in addition to preserving a perfect fluid form of the Brown- York stress energy tensor 
associated to the slice, the timelike Killing vector and the homogeneity in the x % direction 
8|. One transformation is a translation of the radial coordinate, plus a re-scaling of t 

r^r-r h , ^(1- r h /r c )- 1/2 t , (29) 

which moves the horizon from r = to an r = < r c . The second is a boost in the x l 
directions. The resulting metric can be written in the following manifestly boost covariant 
form 

ds 2 = g A Bdx A dx B = -(1 + p 2 (r - r c ))u fl u u dx^dx u - Ipu^dx^dr + P^dx^dx" . (30) 
In this line element we replaced with the relativistic pressure p using the formula 

P= t= ( 31 ) 

V r c ~ r h 

The coordinates x M = (t, x l ). The fluid (boost) velocity is defined as u M = 7(1, v l ), where 
7 = (rl - v 2 y 1/2 and P^ = 7^ + u^u v . 

One can now investigate the hydrodynamic system dual to the above metric. To do 
that, we need to consider the dynamics of the metric perturbations within a hydrodynamic 



limit. In the literature 



BE 



8], the first method developed was to work in terms of the non- 
relativistic variables v 1 and rh and allow these to be functions of space and time: v *(£, x l ) and 
rh{t, x l ), while r c remains fixed. The metric is no longer flat and no longer a solution of the 
vacuum Einstein equation. One then introduces a particular non-relativistic hydrodynamical 
expansion, first proposed in 0, [loj . In terms of a small parameter e, 

v i ~ ev\ex\ e 2 t) P ~ e 2 P(ex\ e 2 t) , (32) 
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where the non-relativistic pressure P(t,x l ) is defined as a small perturbation of the horizon 
radius 

r h = + 2P + O(e 4 ) . (33) 

Expanding the metric (13"U1) out to 0(e 2 ) in this manner yields the non-relativistic metric 
originally found by |z] 

ds 2 = —rdt 2 + 2dtdr + dx;dx % 



2(1 — —] Vidx l dt - —dx l dr 
r r J r r 



r 

+ I 1-- 

rv 



2 



(v 2 + 2P)dt 2 + ^dx l dx J +(— + —] dtdr . (34) 
r c J V r c r c / 

In GR, the momentum constraint equations on the surface 5* c are equivalent to the divergence 

of the Brown- York stress tensor 

R, A N A = d v T*J , (35) 

and constrain the fluid variables. At second and third order in e, momentum constraint 
equations are 

<V = r^» + « 3 )=fl. (36) 

At second order, this is equivalent to the incompressibility condition diV 1 = 0. At third 
order one finds the Navier-Stokes equations with a particular kinematic viscosity 

dtVi + v j djVi + diP - r c d 2 Vi = . (37) 

Imposing the incompressiblity condition implies (1341) is a solution to the Einstein equations 
up to 0(e 3 ). In [sj the higher order corrections to (1341 and the corresponding corrections to 
the incompressibility and Navier-Stokes equations were determined up to 0(e 6 ). 

Compere, et. al also noted that there should be a fully relativistic description of the fluid 
dual to Rindler spacetime. In this case the fluid variables the fluid velocity and pressure, 
which are allowed to depend on 2^: u^{x^) and p(x M ). They constructed the general form 
of a viscous fluid stress tensor to second order in derivatives with zero equilibrium energy 
density. The resulting stress tensor is 

= pu^u u + pP^ - 2t?/C^ 
+ aJCftxv + c 2 /C^| Ak) + csO/Oa, + c A P^P°D x D a hip 
+ c 5 IC^Dlnp + c fi D^lnpDj\np , (38) 



v 



2P\ 
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where D = = P»d u , = P^P^d (x u a) , and = P^P^d [x u a] . The coefficient 77 

at first order is the usual shear viscosity. Note the absence of a bulk viscosity term in ( 138]) 
at the same order. This is due to the fact that at viscous order we can impose the ideal 
order equation d^u^ = 0, which follows from p = 0. On the other hand, there can be viscous 
corrections to the energy density p, which can be parameterized very generally as 

p = bxK^Kr + b^^QT + b 3 D \np D hip + b 4 D 2 hip + b 5 D^ \npD ±fl \np . (39) 

The Cj, i = 1..6, and bj, j = 1..5, are the possible new transport coefficients. Note that 
beyond equilibrium, the pressure and fluid velocity are no longer uniquely defined. This 
ambiguity is usually fixed by a choice of "frame". Usually one works with the so-called 



Landau frame 



201 ]. which is constructed so that the viscous fluid velocity is defined as the 



velocity of energy transport T^u a = and p does not receive viscous corrections. However, 
as we will see explicitly later on, this choice obviously turns out to be inconsistent in the 
Rindler fluid case at second viscous order. Instead one requires 

T$P»u° = , (40) 

and the pressure receives no corrections. 

If one makes a non-relativistic expansion of (!38|) in terms of v % and P and matches to the 
Brown- York stress tensor at 0(e 4 ) and 0(e 5 ) computed using the non-relativistic solution 



found in [8 1, the values of the following transport coefficients can be read off 

77 = 1, bi = -2y/F c , b 2 = 0, ci = -2^/tc, c 3 = -40^, c 2 = c 4 = -4^ . (41) 

However, to fix the remaining second order transport coefficients one has to work to even 
higher orders in the non-relativistic e expansion. Instead we will take a more direct approach 
working with a relativistic hydrodynamic expansion. 

B. The relativistic fluid metric and stress tensor 

In this section we consider perturbations to the metric ( l30l) . treating u^(x^) and p(x^) but 
eaving r c fixed. This follows the standard approach used in the fluid-gravity correspondence 
l|. Now the metric is 

ds 2 = gf B dx A dx B = -^(x^) Ull (x^u u (x^dx^dx u - 2p{x tl )u tl dx' x dr + P^x^dx^dx" , (42) 
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where we have defined $ = l+p 2 (r — r c ) for convenience. This metric is no longer a solution 
to the vacuum Einstein equations, but one can work order by order in an expansion in the 
derivatives of m m and p with respect to d^. This is equivalent physically to an expansion 
in a small, dimensionless Knudsen number A = J2 )^ £ , where £ m f p is the mean free path of 
the underlying system and L the characteristic scale of the perturbations to this system. 
With this solution in hand, we can compute the Brown- York stress tensor order by order. 
This stress tensor should have the form of (1381) . which will allow us to read off the transport 
coefficients in a direct way. 

As a first step to illustrate how this works, we compute the Brown- York stress tensor for 
the metric (j42p at r = r c . This is a solution at zeroth order i.e. Rab = + 0(A). The 
space-like unit normal to this surface is 

Therefore, we need to find the inverse metric (this is also needed in the our calculations later 
on). Using the formula g AC Qcb = Sg, we find 

g" = p- 2 $; g r » = p~ V ; g V» = pV» . ( 44 ) 

Thus, we find n r = p~ x and = . Using 

Kp, = - {n A V A ln» + l^d v n x + jx^n*) , (45) 

(jl]) gives 

T^dafds? = pP llu dx^dx u , (46) 

which as expected is the ideal part of fl38l) with p = 0. 

The strategy for solving the field equations is as follows. One introduces a first order 
correction to the metric Sg^\ 

g = g® + 5gV> . (47) 

The corrected metric at first order induces a 8R^ B at the same order (necessarily involving 
only radial derivatives). We want to solve for the metric 5g^ so that 

5RU + R% = , (48) 
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where R AB comes from the zeroth order metric. This method can be generalized to solve 
for the metrics at higher order in A. If we have a solution to (n — 1) order g^""~ l \ then one 
introduces a correction 5g^ so that 

5Rf B + Rf B = . (49) 
The first step is to compute the general form of SR AB . The Christoffel symbols are 

^ a bc = \g$ (Vb^S + VcSg%l - V D 6g&) , (50) 

where Va is covariant derivative with respect to the background Rindler metric. Hence, we 
need the Christoffel symbols for the Rindler metric. These have the form 

fV = ^$?VV; ?V = ^pu^, f^A = ^pu^u u u x (51) 

with the rest being zero. 

In our solution, we will choose the gauge so that at all orders 

g rr = 0; g ril = -pu^ . (52) 

This implies that 5g$ = and Sgffl = 0. Expanding out ( 1501) gives the following results: 



ST r rr = (53) 

<&V = -p-^driSg^) (54) 

ST\ U = -\p- 2 ®d r {8g$) - l -u,u u u\°5g { £ (55) 

<$T" rr = (56) 

5T\ U = l -P^d r {5g^) (57) 

5V\ X = l - (-p^u^gS) - pP^u v u x u°8g&) (58) 

Now we use the formula 

SRaI = -VaST c cb + V C ST C AB . (59) 
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The final result is 



5R in) = J-^SgS) (60) 
5R^ = IpuMP^SgM) + \v- l dl{uH$) (61) 

- Ui,u v d r {u x u°5g { £) + U>u,u u d r {P x °5g^). (62) 



Notice that these satisfy 



5R ( :y = . (63) 



Using fH9|) we can now obtain the general solution for 5g^ B . This solution needs to be 
consistent with the following boundary conditions: (i) no singularity appearing at r = 
and (ii) induced metric remains flat on r = r c . The second implies that all the n > 1 order 
corrections must vanish at r = r c . Projecting into components normal and transverse to 
we find 

P^SgS = V f ^dr' jf P^Mdr" , (64) 

c 'p 1 



u X P°M£ = (1/2)(1 - r/r c )V^(x) - 2p f' dr' f° dr"P x R^ , (65) 

J r J r' 

u x u a 5g { £ = (1 - r/r c )A {n \x) + p J° dr 1 J*" dr" (pP Xo R^J - p' 1 ®^ - 2fl£V) . 

(66) 

where (V^u^ = 0) and are free, undetermined functions at this stage. 



1. The first viscous order 

With this general solution for any n, we now consider the first order solution, which 
requires the R^b computed from the zeroth order equilibrium Rindler metric. To start, we 
need the Christoffel symbols to first order in A. These can be read off from the second order 
connection presented in Appendix A. 
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Using these results, we find ultimately for Ricci tensor components 

Bg} = 
R§ = 

^JS = d(vP u ») + D P u v u v + ^p{d x u x )u^u u + pu {fl a u) , (67) 

where the operator D = u x d\ and a M = u u d v u^. 

The equation R^An A = is the momentum constraint. However, as we saw earlier ( 1631) 
the SR^a piece satisfies this condition automatically. Therefore we must have u u R^ v = 0. 
Projecting along -u M and orthogonal to with the projector, one finds the relativistic ideal 
hydro equations 

<9X = > (68) 
a M + i*0 A lnp = O. (69) 

So, as expected, the momentum constraints are the relativistic Navier-Stokes equations. 
Inserting ( 157|) into (16411661) yields the solution 

^Pu^l = (1/2)(1 - r/r c )K {1) (^) , (70) 
PtPZSgg = , (71) 
uWSgft = (1 - r/r c )A (1) (a^) . (72) 

To determine these functions, we need to impose the "frame" conditions on the 1st order, 
viscous part of the Brown- York stress tensor. Using we find the first order part of the 
stress tensor is 

T$ = ({tcPVA^ + 2p- l Dp) j, v + (r c py\A^u,u v - u ifl V^) 

- 2d^u v) - 2p~ 1 u^d v) p, (73) 

where we have imposed d^v? = 0. We now require the stress tensor satisfies the Landau-like 
condition (j4"0l) . This yields an equation for V^: 

V® = -2r c P;d v p + 2r c pa M . (74) 

The second condition we demand is for p to be the pressure at all viscous orders. Thus the 
term proportional to 7^ must vanish. This implies 

A {1) = -2r c pD(\np) . (75) 
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Feeding these results back into the stress tensor, we find simply 

2$ = -2/C„ . (76) 

So, as we expected, the shear viscosity 77 = 1. The complete solution for the metric to first 
order is (imposing the ideal hydrodynamics equations) 

ds 2 = -(1 + p 2 (r - r c ))u^u u dx^dx v - 2pu p dx p dr + P pv dx p dx u 

+ 2p(r — r c )D(\np)u p u v dx p 'dx u — Ap(r — r c )u^P^d\ \npdx p dx u . (77) 

2. The second viscous order 

In order to solve to second order in A, we need to find pS^ b . The first step is to compute 
the connections out to second order. Then it's a matter of grinding though the calculations 
of the Ricci tensor, determining the solution, and then computing the Brown- York stress 
tensor. This will fix all of the second order q transport coefficients. The connections and the 
Ricci tensor that come from metric (ITT)) are complicated, so we present them in Appendix 
A. 

We find that at second order the momentum constraint equations R^ati a = projected 
once on u p and once on P£ are: 

- p~ 1 d a u x d x u' T - p- l P pa d p u u d a u u = (78) 
CV + P x d x lnp - p^P^d^ + p- l u p P va d v u x d a u x + 2p- x P vu d vUil dMp = 0. (79) 

Solving the equations (I64H66P for Sg^J we get: 

u"P}5g$ = l -p 2 {r - r c ) 2 {-2P^u p d p lnp + 2P Xp d x u v d p \np - P pa d p d a u v + u u P^d a u p d p u p ) 

+ 5(l-f)VjV) (80) 

PiP^gfl = (r- r c )P«P? {-dMpdfMp + 2d a d^vp + 2d {a u p) u p d p \np - 2u^d a d {a u^ 

3 1 

+2u a d a u^d a) \np - a a a - -d a u x d p u x + -P ap d a u [S d p u a + d a u {a d 0) u 

+ \p\r - r c ) 2 P«P? (-2d a u {a d^u a + P° x 8 a upd x u a + d a u%u x ) (81) 

u p u y 5gfi> = ^p 2 {r - r c ) 2 (d a u a d a u a + P° x d a u v d x u v + AP^dJnpd^np) 

+ ~p 4 (r - r c ) 3 (-d a u a d a u a + P aX d a upd xU P) + (1 - -)A^ 2 \x p ). (82) 
4 r c 
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The corresponding Brown- York stress energy tensor is: 

Tpv = 7a* (P + P" 1 (^"^ (2) + ^n A a A n CT + P^d^d^)) + pu P u u 
+ (-29( jU )W I/ ) + 2w (M <9 i ,)hrp + 2u p u u u x d x lnp) 

- p^P^pP {-Ad a Yn.pdp\n.p + 2(9 a (9arnp + 2d( a u l 3)U p d p \n.p - 2u a d a d( OL Ui3) 

3 1 \ 

--d a u x d p u x + -P ap d a upd p u a + d^u^d^u J . (83) 

As before, we impose the frame condition that the pressure be unchanged from its equilibrium 
value. This condition eliminates the derivative parts of first term proportional to 7^ in fl83|) 
and fixes to a non-zero value. Note however that the stress tensor now has a term 
proportional to u p u v (in the third line above). Hence at this point we see the traditional 
Landau frame T pv u p = will be inconsistent. To fix we instead require ( I40j) . 
With these values fixed, the stress tensor can be put in a more conventional form: 

T pu = pP^ v - 2/C M „ - 2p- 1 u p u v /C a/3 /C a/3 - 2p~ x KJCl - Ap' 1 ^^ - Ap^Q^ 
- Ap- x P«pPd a d p hip - Ap^K^Dlnp + Ap' 1 D^lnpD^hip (84) 

Due to the u^Uy term, one can see that the energy density is no longer zero. It gets following 
correction in the second order: 

p = T pv u»u v = -p^d^dxv? - p- l P p °d p u v d a u v = -hc^. (85) 

We can also read of from the stress tensor the transport coefficients. In the notation of js], 
we get: 

ci = -2p _1 , c 2 = c 3 = c 4 = c 5 = -4p _1 , c 6 = 4p -1 . (86) 

which is in agreement with [8J, who found the first four transport coefficients and the energy 
density. 
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(2) 

Finally for the metric solution og)/J we get: 

6g$ = u»u v Qp 2 (r - r c ) 2 {lC a ^ + 2P° x dJnpd x \np) 

+ \v\r - r c ) 3 (O^) + 2(r - r c ) (/C^/C^)) 
+ Qp 2 (r - r c ) 2 (4ft p CT «9 CT hrp + P Xa d x d a u p ) 

-{r - r c )(-P x °d x d a u p + 2K;dJnp - 2Q;dJnp) 

+ (r - r c ) + 4/Cf M fi| pH + 4n w n" v + AP^d a d p \np + 

-4/^lnp/^lnp) 

+ p 2 (r - r c f (0 MP 0/) . (87) 
For reference, the inverse metric of (l87|) is presented in Appendix B. 

IV. THE ENTROPY CURRENT 



In 14] it was shown that in the fluid-gravity correspondence the entropy current of the 
dual fluid on the boundary can be mapped into the area current of the black hole event 
horizon. Thus the second law of thermodynamics is equivalent on a geometrical level to 
Hawking's area theorem. Here we will follow this general prescription to calculate the 
entropy current for the Rindler fluid to second order in the gradient expansion. Given the 
exotic properties of the Rindler fluid, it is clearly of interest to determine whether it behaves 
consistently with the second law. 

First, since metric solution is no longer stationary, the event horizon is dynamical and its 
location varies in time and space. In order to find r/ l (x M ) we will need to solve the following 
equation in the derivative expansion 

g AB d A (r - r h {x»))d B {r - r h {x»)) = . (88) 

Using our previous results for the metric, it is straightforward to show at second order 



r h = r c - 1 + -u»d,lnp - ^JC^JC^ - 

p l p 6 2p 4 2p^ 



8 1 4 

- — DhipDhip + —D^lnpD^kip + —D(D\np) . (89) 
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We can define a co-dimension 2 hyper-surface by two null normals to the hyper- surface. The 
ingoing null geodesies n A and the outgoing null geodesies I are: 



£» = A", £ r = B, n r = -1, = . (90) 

The unknown functions N 1 and B can be found from the following relations: 

£ A £ A = 0, n A n A = 0, £ A n A = -1 . (91) 

From these conditions we see that there is another freedom in determining £ M in any order 
except from the zeroth order. Therefore, we impose the requirement that the vector £ A will 
be the normal vector to foliations of hyper-surfaces that do not intersect with each other. 
This requirement is called the Frobenius condition and is expressed by the equation: 

v A dv, v A = 9ab£ B . (92) 

This gives us, on the event horizon (putting r = r/ l (x M )) , the same null normal vector that 
we get by calculating the normal to the event horizon directly from the equation that defines 
the normal to the event horizon: £ A = g AB ds(r — r^). 
To second order, the vector £ A on the horizon is: 

2 6 2 2 
£ r = Fd.rh = —u^djnp -M /i <9 M hTpu p cUnp -P^ p djnpdp\iap + —u^u p dud p hip (93) 

pi pA p4 

F = -u p - —^ p P Xa d x d a u p - —^ x d x \np + ^P> lu d u \npu p d p \np - ^P^u p d u d p \np . 
p 2p A p 6 p 6 p d 

(94) 

In order to compute the entropy current we will employ horizon expansion 8(e) along the 
horizon generator t A 

0(i) = (g AB + i A n B + t B n A ) V A £ B = l -d, (V^ (0> + Vgp^ + V9~pt {2) ») ■ (95) 
We can identify the entropy current as the term in the brackets up to an overall factor of 



1/4G (which in our units of 16%G = 1 is An) [21]. 



In order to compute the entropy current we therefore need two ingredients: The square root 
of the metric determinant, and the null generator £ p . We explained how to get the latter. 
The former can be derived by computing the expansion of the null normal £ A , where £ r = £ r 
and £ u = £^ ^, then we will get only the first term in the brackets of (f95l) and we can identify 
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immediately the square root of the determinant of the metric. 

Combining all the ingredients we get the following result for entropy current: 

S» = ^ (l - j 2 (lC a ^ - hl a ^ + 2P a ^d a d \n P + 2/CDlnp - 2P^d a \npd^ ) . 

(96) 

Taking the divergence and imposing the Navier-Stokes equations gives 

dpS* = ^ (Kap + l - {-ZKapu^hxp + 2d a \npdp\np - 2d a d p \np - 2/C£/C M/3 - 2fi/ft^ 

(97) 

which is clearly non- negative, just as expected from the area increase theorem applied to 
the Rindler horizon. 



V. DISCUSSION 



In this paper we extended the fluid/gravity correspondence to relativistic Lorentzian ge- 
ometries. The fluid is defined on a codimension one timelike hypersurface, embedded in 
the bulk space-time and equipped with a flat intrinsic metric. We considered two types of 
geometries that are solutions to the vacuum Einstein equations: the Rindler metric and the 
Taub plane symmetric vacuum. They are found correspond to dual perfect fluids with van- 
ishing and negative energy densities, respectively. The Rindler geometry is characterized by 
a causal horizon and one can systematically construct its hydrodynamic derivative expan- 
sion. The Taub geometry, however, has a timelike naked singularity, indicating pathological 
behavior. 

The dual fluid in the Rindler case has a vanishing energy density as its equation of state. 
While such an equation of state is rather unusual, the hydrodynamic expansion up to second 
order exhibits no pathologies and the gradient of the entropy current is still positive. We 
determined the transport coefficients and verified their consistency with those obtained in 
the non-relativistic Rindler hydrodynamics. 

There are various directions to proceed from here. One may wish to consider the higher 
curvature corrections the fluid/Rindler relation. In the AdS/CFT examples this provides 
information on the holographic relation to field theory deformations as well as at subleading 
orders in the field theory strong coupling expansion. Also here, such an analysis is likely 
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to add valuable information on the correspondence. In particular, now the entropy current 



is no longer the area current 22| and one needs to define it appropriately and verify the 
positivity of its gradient. One may also wish to consider the field theory dual to Rindler 
geometry and study its properties. There is at least one candidate non-local dual field theory 
8|, which gives the correct equation of state at the ideal fluid order. One can try to study 
its thermal properties and its hydrodynamic limit. This can shed light on holography in 
asymptotically flat spaces, where the dual field theory is expected to be a non-local one. 



Note added: Sections III and IV contain some overlap with [23j which is posted simul- 
taneously with this paper on the ArXiv. 

Appendix A: The Connections and Ricci Tensor 

The connections up to 2nd order from the first order metric (I77p are: 

r V = ^( 3 <Vm9 + pu p + u p u x d x \np + u a d a u p ) 

+ (r - r c ) (d p u x d x \np + u p P xp d p \npd x \np - P xp d p u p d x \np) (Al) 

I i 

TV = + U(yd v) h\p + -pu^Uv + u p u v u x d x \np) 

p z 

+ p(r - r c )(3u {p d u) \np + ^pu p u u + ^u x d x (u p u u ) + u p u u u x d x \np) 

+ p 2 (r - r c ) 2 (2u( M d u )U p d p \np - P Xp d x (u p u u )d p lnp) 

+ (r - r c ) [2d l Mpd y hip + 2d p d u \iap + 2u^d p )U p d p \np + 2d( p u u )U p d p hip 

— u p u„u x d x \npu p d p \np + u x d x (u p u u )u p d p \np + u p u u u x d x u p d p \np 
+u p u u u p u x d p d x lnp + 2u x d x u^d u )lnp + 2u( p u x d x d u )\np) (A2) 
r V = \v - ^P^lnp + P^cUj (A3) 

Y v pa = - (pu u u p u a - 2u v d( p u a ) + 2u v u {<T d^\np + 2u v u p u (7 u x d x \np) 
+ \v\r - r c ) {-2u {a d A u v + P vX d x {u c u„)) 
+ p(r - r c )P up d p lnp (2d {p u a) - 2u^d a) \np - 2u p u a u x d x \np) 
+ p{r - r c )P uX (2u( a d^u x u p d p \np - 2d^u x P p ) d p \np - 2d (ji (pu a )P£d p lnp) 

-d x {u a u p u p d p \np - 2u ip P p ) d p \np)) (A4) 
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The Ricci tensor calculated from the 1st order metric (IT?]) up to 2nd order is: 

RW = ^p(2u p u p d p \np - 2P» p d v u p d p \np + u p P uX d p \npdMp - d v d^u v - u p d u u u u x d x \np 
- u p P uX d u d x \np + cU*VcU M + P v °d v d a u^) 

+ \p\r - r^-u^dvu" + u^dvu'dxu,) (A5) 

1 3 
R< $ = --dfMpd u lnp + 2d p d u lnp + 2u( p d u )U p d p lnp + d^u v )U p d p \ap + -u x d(u p u u )u p d p lnp 

+ 2u p u u u x d x u p d p lnp + 2u p u u u p u x d x d p \xip + 2u^u x d x d u )lnp - d a u a d^u v) - u u d a d^u u ) 
+ d a u a U( p d u )lnp + U( p u a d a d u )\np + u a d a u^d u )lnp + u p u u d a u a u x d x lnp - d^hip 
- U( p d u )\npu x d x \iap - -u Pl u u u p d p lnpu x d x lnp - ^u p d p u p u x d x u u - ^d p u a d u u a 
+ ^P ap d a u v d pUfl - ^P° p d p (u p u„)djnp + ^UvP^dMpdxlnp 

+ p 2 {r - r c ) (-d a U( u d p )U a - u^ u d p )d a u a + ^d a u a u x d x (u p u u ) + u {p d u) u° "<9 CT lnp 

1 1 

- -P ap d p (u p u u )dMp + u p u v u x d x u a dMp + -d p u x d u u x - P ap d p u u d a u^ 

+ \p\r - r c ) 2 ( u Pi u u 3 x u <T 3 cr u x + Ufi u v P ap d a u x d p u x ) . (A6) 



Appendix B: Inverse metric 



The inverse metric of (1871) to the 2nd order is: 



1 2 

9 rr = (l + P 2 ( r - r c)) (r — r c )u x d x \np 

p z p 

- ~(r - r c ) 2 {K a pK^ - 6P pX d x \npd p \np) - X -p\r - r c ) 3 Q a ^ 

-^{r-r c )K aP K^ (Bl) 

g T v = - u p - 2(r - r c )P px d x \np + -p(r - r c ) 2 UQ p(7 dJnp + P pp P Xa d x d a u p ) 
p 2 

--{r- r c ) {-P pp P Xu d x d a u p + 2K pa dMp ~ 2Q pa dJnp) (B2) 
gT = - (2K, P 1C P + 4JC p ^n^ + AQ p p Q pu + AP pQ P u ^d a d lnp + AIC pu Dlnp 

-AD ±p lnpD ±u lnp) (r - r c ) - p 2 {r - r c ) 2 Q p p Q pv (B3) 
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